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SUMMARY 


Sonic boom propagation is not steady, even in a nonturbulent atmosphere. The shock 
overpressure A p is affected by geometrical spreading, stratification of the atmosphere, and even the 
N shape of the waveform. Nevertheless, for purposes of predicting shock profile and rise time, it has 
commonly been assumed that the shock is in steady state. For example, molecular relaxation, which 
is a major factor controlling sonic boom rise time, is strongly dependent on relative humidity. 
Because humidity varies with altitude, rise time varies as the sonic boom propagates downward. 

The question is whether rise time depends only on local conditions or is also affected by the history 
of humidity variation along the propagation path. Kang [1, Chap. 7.2] argues that shocks respond 
to change in humidity quickly enough that they are in effect always in steady state. In other words 
only local conditions are important. Robinson [2, Chap 5.2] however disagrees with this hypothesis. 
Raspet et. al. [3] found that perturbed 100 Pa shocks (step waveform) require propagation 
distances of order 1 km for the rise time to return to within 10% of its steady shock value. 

The purpose of our investigation is to determine the effect of unsteadiness (not associated with 
turbulence) on rise time. The unsteadiness considered here is due to (1) geometrical spreading, (2) 
stratification, which includes variation in density, temperature, and relative humidity, and (3) N 
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shaped waveform. A very general Burgers equation, which includes all these effects, is the 
propagation model for our study. The equation is solved by a new computational algorithm in 
which all the calculations are done in the time domain. 

The present paper is a progress report in which some of the factors contributing to unsteadiness 
are studied, namely geometrical spreading and variation in relative humidity. The work of Pierce 
and Kang [4], which motivated our study, is first reviewed. We proceed with a discussion of the 
Burgers equation model and the algorithm for solving the equation. Some comparison tests to 
establish the validity of the algorithm are presented. The algorithm is then used to determine the 
distance required for a steady-state shock, on encountering an abrupt change in relative humidity, 
to reach a new steady state based on the new humidity. It is found that the transition distance for 
plane shocks of amplitude 70 Pa is about 4 km when the change in relative humidity is 10%. Shocks 
of amplitude 140 Pa require less distance. The effect of spherical and cylindrical spreading is also 
considered. We demonstrate that a spreading shock wave never reaches steady state and that its 
rise time will be less than the equivalent steady state shock. Finally we show that an N wave has a 
slightly short rise time than a step shock of the same amplitude. 
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1. Change in relative humidity 

2. Geometrical spreading 

3. Waveform 


REVIEW OF THE PIERCE-KANG PREDICTION OF SONIC BOOM RISE TIME 


Pierce and Kang [4] calculated sonic boom rise time by solving a model propagation equation 
called Burgers’ equation [5]. An important contribution was to augment the classical Burgers 
equation with terms that describe molecular relaxation due to both nitrogen and oxygen. By 
assuming that the sonic boom shock near the ground is in steady state, Pierce and Kang simplified 
their Burgers equation and solved it numerically. Rise time predictions obtained from their solution 
were then compared with rise time data from a large number of sonic boom measurements taken at 
Edwards Air Force Base in 1987 [6]. In their initial comparison almost all the measured rise times 
lay well above the predicted curve [7], as shown in the left-hand plot of Fig. 1. The discrepancy was 
attributed to atmospheric turbulence. Later, however, when they made a correction for pressure 
doubling at the ground, they found their prediction to fall in the middle of the data [8]; see the 
right-hand plot in Fig. 1. Although the data would seem to corroborate their prediction, the role of 
turbulence casts doubt on this conclusion. Model experiments done at our laboratory [9, 10] 
indicate that turbulence almost always increases rise time, rarely decreases it. If turbulence does 
have this effect, then molecular relaxation should be expected to provide a lower bound for the 
data, not approximate a mean for the data. Why then does the Pierce-Kang prediction not serve as 
a lower bound? 



Steady state shock overpressure (Pa) 



Figure 1: The left-hand plot shows Kang and Pierce’s initial prediction [7]. The right-hand plot 
shows the corrected prediction [8]. 
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Our hypothesis to answer this question is based on results found in Refs. [2, 3]. We suspect that 
the shock wave at the head of a sonic boom does not respond quickly enough to variation in 
atmospheric conditions (and to other changes that affect the profile, such as geometrical spreading 
and even wave shape), to justify the steady-state assumption. If our hypothesis is correct, then to 
improve on the Pierce-Kang prediction requires that more than local conditions be taken into 
account. Past history along the propagation path must be significant. Figure 2 shows profiles of 
temperature and relative humidity measured during a sonic boom experiment [1, pp 157-161]. It is 
seen that conditions can change rapidly, particularly during the lower part of the propagation path. 
To determine whether the sonic boom profile can respond quickly to changes of this order, we have 
calculated the effect on rise time of an abrupt change in atmospheric conditions. Thus far we have 
concentrated on changes in relative humidity. For purposes of this paper, temperature and pressure 
are fixed at their ground level values. 




Figure 2: Atmospheric conditions measured during the sonic boom experiment. 
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THE BURGERS EQUATION 


The “classical” Burgers equation is the standard model equation for plane finite-amplitude waves 
in a thermoviscous medium: 


d 2 p 

dx 2 p 0 c$ dt' TV dt' 2 ' W 

Here p is acoustic pressure, t time, t' = t ~ x/co retarded time, x distance, Co small-signal sound 
speed, po ambient density, (3 coefficient of nonlinearity, and S T v the thermoviscous loss coefficient. 
Pierce [11] added terms to account for relaxation processes. Each relaxation process v is 
characterized by a relaxation time r u and a change in small-signal sound speed (A c)„ due to the 
relaxation. In operator notation Pierce’s “augmented Burgers equation” may be written 
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Equation 2 is still for plane waves. If geometrical spreading is included, the equation becomes 
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where the spreading factor a is 0 for plane waves, \ for cylindrical waves, and 1 for spherical waves 
This is the equation we have solved numerically to obtain the results reported in this paper. 


Burgers’ equation may be further generalized to include effects of (1) stratification and (2) 
diffraction. Stratification may be included by scaling the dependent variable p and stretching the 
independent range variable x [12]. The Burgers equation for this case is not considered in this 
report but will be taken up in the future. To include diffraction effects, one must use the KZK 
equation, which is a multi-dimensional form of Burgers’ equation; see, for example, Refs. [13, 14]. 
As a spinoff from the present work, relaxation effects have been included in a computer code that 
solves the KZK equation, but no formal report of the results has yet been given. 
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NUMERICAL PROPAGATION 


Solutions of the generalized Burgers equation that are not in steady state involve solving a 
partial differential equation. Except for a few rare cases the solution can only be obtained 
numerically and it is common to use some sort of marching scheme. A time waveform is digitized 
with M samples and then small steps are taken in the propagation direction. At each step 
absorption and nonlinearity are solved in series. It is popular to do the absorption effects in the 
frequency domain as this requires M complex multiplications. However in the frequency domain the 
nonlinear term involves a convolution — which requires of the order of M 2 operations. To speed up 
the code the nonlinear distortion can be applied in the time domain as it requires only order M 
operations. To go between the time and frequency domain one can use the fast fourier transform 
which requires order M\ogM operations. Algorithms like the Pestorius [15] code flip-flop between 
the time and frequency domains at each step to take advantage of computing absorption in the 
frequency domain and nonlinear distortion in the time domain, the price being the use of the FFT. 


X = X4-AY 



Nonlinearity 


Figure 3: Pestorius type approach to solving the Burgers equation. 
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TIME DOMAIN ALGORITHM 


It would be nice to stay in one domain but without having to pay the price of M 2 operations for 
a convolution. Lee and Hamilton [13, 14] have developed a method of computing the absorption in 
the time domain. They approximate the absorption with a finite difference equation. This yields a 
tn diagonal matrix system which can be solved in order M operations. The code they developed was 
used to solve the KZK equation, which is a generalization of the Burgers equation to include 
diffraction effects. We have extended the code to account for molecular relaxation and spreading 
effects. The work presented here does not, however, include diffraction effects. 

Apart from its numerical advantage the fully time domain algorithm has the nice property that it 
can propagate pulses. Because the FFT isn’t used, the endpoints of the waveform need not match 
to make a periodic waveform. Therefore step shocks and N waves are easily dealt with. The 
algorithm is particularly suited to finding a steady state solution. Raspet et. al [3] had to use a 
square pulse to find the steady state behavior of a shock. Square pulses have a limited propagation 
range before they turn into sawtooth waves. In the time domain code a pure shock can be 
propagated with out difficulty. 
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Figure 4: Time domain approach to solving the Burgers equation. 
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STEADY STATE: NONLINEARITY AND THERMOVISCOUS LOSSES 

A number of cases were run to test the validity of the code. The first was to obtain the 
steady-state solution of the classical Burgers equation for a thermoviscous fluid. The known 
analytical solution for the steady shock is the hyperbolic tangent function. Figure 5 shows how a 
shock front is propagated with the time domain code; u is the distance variable. The first figure 
shows the initial profile, chosen because it looked interesting. The other figures show how the profile 
develops. The final figure, a = 2, shows that the numerical result agrees very well with the 
analytical solution. 



-200 0 200 -200 0 200 
t' t’ 

o = 1 o = 2 



Figure 5: Propagation of a shock in a thermoviscous medium. 
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STEADY STATE: NONLINEARITY AND ONE RELAXATION PROCESS 


We now determine whether the relaxation part of the code behaves correctly. Polykova et. al. 
[16] obtained the steady-state solution for a finite amplitude wave in a medium with one relaxation 
process but no thermoviscous losses. Their result (denoted PSK in Fig 6) is 


_ i n C 1 + P/po) D 1 
T u (I+P/Po) 0 ' 1 ' 1 ’ 


where 
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Pop 

Figure 6 shows the result from the propagation program in a monorelaxing fluid. For the values 
chosen relaxation was not enough to stop the waveform from becoming multivalued. In the 
analytical result weak shock theory was used to ensure a single valued function. Multivaluedness 
was prevented in the numerical algorithm by including a small amount of thermoviscous losses. 



Figure 6: The top figure shows the analytical result for the steady-state solution in relaxing medium 
with no thermoviscous effects; D ~ 0.5. The middle figure shows the initial and steady-state profiles 
obtained by the time domain code. The bottom figure compares the analytical and numerical steady- 
state profiles. 
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STEADY STATE IN AIR 


The last verification test was against Kang’s result. A plane shock front is sent into a standard 
atmosphere with a relative humidity of 10%. This allows us to compare results with a steady-state 
result in Kang’s thesis. In our calculation the shock started out with a hyperbolic tangent profile 
and was then propagated until the profile no longer changed. Figure 7 compares the two results. 



lime relative ui shock arrival (ms) 



Figure 7: Steady-state solution in air; T — 20 °C, Pq — 1 atm, and a relative humidity of 10%. On 
the left is Kang’s profile [1, Fig. 5.8]. On the right is the profile from the time domain code. 
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EFFECTS ON RISE TIME 


We will now use the time domain code to investigate the behavior of the rise time of shocks in 
the air. In all calculations the temperature is 20 °C, and the pressure is 1 atm. 

First we examine how long it takes waveform to recover from a small but abrupt change in 
i via live humidity. A plane wave shock is propagated in air of given relative humidity until it reaches 
steady state. This steady-state waveform is then used as the input waveform for an atmosphere 
with another relative humidity. 


Second we investigate the effect of spreading on the rise time of a shock front. Shocks that are in 
steady state are propagated as spreading waves; because the amplitude of the shock decreases the 
rise time is expected to increase. We present results for the combined effects of spreading and a 

change m relative humidity. Finally we investigate the difference between the rise time of an N wave 
and a step shock. 


TRANSITION DISTANCE FOR 70 Pa SHOCKS 


We use the term transition distance to describe how far a shock needs to travel to go from one 
steady-state profile to another. A somewhat similar term “healing distance” is commonly used in 
literature related to turbulence for the distance a perturbed shock needs to return to its original 
state [3]. In this case we shall look at transition distances due to a change in relative humidity. 

Figure 8 shows rise time as a function of propagation distance for a plane step shock of 
amplitude 70 Pa which starts in a medium of 20% relative humidity. The relative humidity of the 
second atmosphere is 10%, 20% or 30%. The results show the transition distance to be at least 
5 km. The plot in Fig. 9 shows the rise time for a shock initially in an atmosphere of 50% relative 
humidity. Transition distances are greater than 2 km. 

The initial fluctuations in the rise time are due to rather gross changes in the profile. The 
changes are such that the 10% to 90% definition is not a very suitable measure of rise time. Similar 
fluctuations were observed by Raspet et. al. [3], 
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Figure 9: Change in rise time for a waveform leaving a medium of 50% relative humidity. 




TRANSITION DISTANCE FOR 140 Pa SHOCKS 


Figure 10 shows the behavior of a shock wave of amplitude 140 Pa which starts off in air of 20% 
relative humidity. The transition distance needed when the new relative humidity is 30% is about 2 
km. However the transition distance is at least 6 km when the new relative humidity is 10%. 



Figure 10: Change in rise time for a 140 Pa shock leaving a medium of 20% relative humidity. 
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SPREADING 


The amplitude of a spherically spreading shock should decrease as 

A p - ~Ap 0 , (4) 

x 

and for cylindrical spreading, 

Ap 0 . 

The steady state rise time found from the analytical solution of the classical Burgers equation is 

Ar = ln(9)%f^ (5) 

As the amplitude of the waveform decreases the rise time increases because the nonlinear steepening 
effects are weaker. However it is not clear that a spreading waveform will be in steady state. This 
would require the absorption mechanism to respond immediately to the spreading. NaugoFnykh [18] 
argued that a spreading shock in a thermoviscous medium should have a rise time that is shorter 
than the steady-state value because the absorption mechanism can’t work fast enough. 

If a spreading shock remains in steady state then, from Eq. (5) the rise time should vary 
inversely as the pressure jump. Since for spherically spreading waves the pressure varies inversely 
with distance, Eq. (4), we have 

At oc x . 

To investigate the validity of this relation we started with the hyperbolic tangent profile appropriate 
for a plane step shock. The shock was then propagated as a spreading wave. Figure 11 shows the 
initial waveform and how the shock diffuses as it loses amplitude. Figure 12 compares the steady 
state prediction of the rise time to the numerically calculated rise time. Note that the steady state 
prediction always overestimates the rise time. Absorption cannot act quickly enough to diffuse the 
profile before more amplitude decrease, due to spreading, occurs. 
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Figure 11: A shock front that starts off in steady-state in a thermoviscous medium is propagated as 
a spherical wave. The upper plot shows the initial waveform. The lower plot shows the waveform at 
20 times the source radius; the steady-state waveform is also shown. 





Spherical Spreading 


Cylindrical Spreading 




Increase Amplitude by 4 




Figure 12: The rise time of a step shock in a thermoviscous medium. The waveform starts off in 
steady-state and is propagated as either a spherically or cylindrically spreading wave. Two different 
initial amplitudes are used. Rise time is normalized to the initial rise time and distance is normalized 
to the source radius. 
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EFFECT OF SPREADING ON TRANSITION DISTANCE 


We now investigate the combined effects of spreading and a change in relative humidity. A source 
is assumed to be 15 km away and the shock front is propagated at an angle corresponding to a sonic 
cone from an aircraft flying at Mach 2. Cylindrical spreading is used. Figure 13 shows the rise time 
curves and compares them to the plane wave curves taken from Fig. 8. In Fig. 14 we compare the 
rise time curves to those of spreading waves in a homogeneous atmosphere, which are the effective 

asymptotes. We see that the waveform is never in steady state but that the asymptotes are reached 
in 6-8 km. 



Figure 13: The rise time of a shock that is cylindrically spreading; the shock starts from a steady-state 
profile for an atmosphere with 20% relative humidity. The dotted lines show the results for plane 
waves as shown in Fig. 8. Note that the steady-state rise time is not achieved. 



Figure 14: The same curves as in Fig. 13, except that the dotted lines show the results for cylindrically 
spreading waves in a homogeneous atmosphere; these are the asymptotic results. 
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EFFECT OF PROFILE ON THE RISE TIME 


Finally we compare rise time for an N wave with that for a steady-state step shock. Figure 15 
shows the curves for a plane N wave that has propagated a long distance in a medium with a 
relative humidity of 50%. The amplitude of the N wave is 48 Pa. We see that a steady-state shock 
of the same amplitude has a longer rise time. However a steady-state shock of amplitude 50.4 Pa 
has the same rise time. It appears that the profile corresponds to that of a higher amplitude shock 
because the amplitude of the N wave decreases with propagation distance. Just as was seen with 
spreading waves, the absorption cannot keep up with the loss in shock amplitude. 



Figure 15: The solid line is an N wave with duration of 125 ms. It has propagated for 30 km. The 
dashed line is a step-shock of the same amplitude. The dot-dash line is a shock of the same rise time 
as the N wave. This shock has an amplitude of 50.4 Pa. 
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CONCLUSION 


Our study indicates that assuming a steady-state profile for the head shock of a sonic boom 
tends to cause an overestimate of the rise time. First, we have found that substantial propagation 
distance is required for a shock to respond to abrupt but small changes in relative humidity. When 
the absorption increases, the actual rise time is always less than the steady-shock value. Second, 
neglect of geometrical spreading results in an overestimate of rise time. Third, modeling the N wave 
as a step shock also raises the rise time, although the effect is small. All of these results are 
consistent with the conclusion that for a decaying shock dissipation is always “struggling to catch 
up with nonlinearity.” 
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